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Abstract. Results from the photon multiplicity measurements using a fine granularity preshower 
photon multiplicity detector (PMD) at CERN SPS are discussed. These include study of pseudo- 
rapidity distributions of photons, scaling of photon multiplicity with number of participating nucle- 
ons, centrality dependence of < pr > of photons, event-by-event fluctuations in photon multiplicity 
and localised charged-neutral fluctuations. Basic features of the PMD to be used in STAR experiment 
at RHIC and in ALICE experiment at LHC are also discussed. 
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1. Introduction 

Measurement of photon multiplicity in relativistic heavy-ion collisions is complimentary 
to the well established methods of charged hadron measurements and it shows a great 
promise in studying the various aspects of the reaction mechanism of phase transition from 
hadronic matter to Quark-Gluon Plasma and dynamics of particle production. In heavy-ion 
collisions it is important to correlate information obtained from various global observables 
(such as charged particle multiplicity, mean transverse momentum (< p? >) and trans- 
verse energy (Et) ) for proper understanding of the physics processes occurring in the 
reaction. Photon multiplicity is an additional global observable. More specifically photon 
multiplicity provides an unique opportunity to study the changes in the relative population 
of the electromagnetic and hadronic components of the multi-particle final state. These 
and several other aspects are discussed in this paper. 

The paper is organised as follows, first we give a brief description of the photon mul- 
tiplicity detectors (PMD) used at SPS experiments. In section 3 we discuss the various 
physics issues addressed by the PMD at SPS, such as the study of pseudo-rapidity distri- 
bution of photons (dNy/drf), scaling of total number of photons with number of partici- 
pating nucleons, measurement of event-by-event mean transverse momentum (< pr >) of 
photons, event-by-event fluctuation in photon multiplicity and search for domains of disori- 
ented chiral condensates (DCC) through the study of localised charged-neutral fluctuations. 
A comparison to the results from charged particles will be made wherever possible. The 
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study of azimuthal anisotropy of photons is covered in an separate article in this volume. 
In section 4 we discuss the PMD to be installed in STAR experiment at RHIC and ALICE 
experiment at LHC. This is followed by a summary in the last section. 

2. Photon multiplicity detectors at CERN SPS 

A fine granularity preshower PMD was implemented in the WA93 experiment at CERN 
SPS [], allowing the first ever measurement of multiplicity, rapidity and azimuthal dis- 
tributions of photons in ultra-relativistic heavy-ion collisions. The basic features of the 
detector are given in Table 1 . The minimum bias distribution of the photon multiplicity as 
measured by the PMD for S + Au at 200 AGeV in WA93 experiment is shown in Figure |l|. 
The distribution has been obtained for the full azimuthal coverage of the PMD and has 
been compared to results obtained from the VENUS event generator. It is observed that 
VENUS under predicts photon multiplicity for central events. Subsequently a similar PMD 
was also installed in the WA98 experiment []. The basic features of the PMD in WA98 ex- 
periment are also given in Table 1 . The minimum bias distribution of photon multiplicity 
as measured by the PMD in WA98 experiment is shown in Figure [l] The experimental re- 
sults for different target ions has been compared to those obtained from the VENUS event 
generator. VENUS is found to under predict photon multiplicity for central collisions, and 
it under predicts more for asymmetric ion collisions (Pb + Nb and Pb + Ni). 



Table 1. PMD in WA93 and WA98 experiments 



Basic features 


WA93 


WA98 


Data taking 


1991 and 1992 


1994-1996 


Beam and target 


S + Au & S + S 


Pb + Pb, Pb + Ni & Pb + Nb 


Energy 


200 A GeV 


158 AGeV 


No. of scintillator pads 


7600 


53000 


Readout 


CCD camera 


CCD camera 


Distance from target 


10.09 m 


21.5 m 


r\ coverage 


2.8 - 5.2 


2.5-4.2 


rj coverage with full </> 


3.3-4.8 


3.2-4.0 


Pt acceptance 


> 20 MeV/c 


> 30 MeV/c 


Photon counting Efficiency (central to peripheral) 


65 - 75% 


68 - 73 % 


Purity of photon sample(central to peripheral) 


70% 


65 - 54 % 



Both WA93 and WA98 experiments include several other detectors for charged parti- 
cle multiplicity and momentum measurements along with trigger detectors for defining 
the centrality of the reaction. The common coverage of the charged particle multiplicity 
detector and the calorimeter (for measuring E^) with PMD was used for event-by-event 
measurement of < px > of photons and search for possible formation of DCC. 

3. Physics issues addressed by PMD at SPS 

In this section we discuss the various physics issues addressed by the PMD at CERN SPS. 
These includes study of dNj/dr/, scaling of photons, event-by-event < px > of photons, 
event-by-event fluctuation in photon multiplicity and formation of DCC. 
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Figure 1. Minimum bias inclusive photon cross sections for S+Au reaction at 
200-AGeV (left plot) and for Pb+Ni, Pb+Nb, and Pb+Pb reactions at 158 AGeV (right 
plot). Solid histograms are the corresponding distributions obtained from the VENUS 
event generator. 



3.1 Pseudo-rapidity distribution of photons 

One of the challenges in relativistic heavy-ion collisions is the large number of particles 
produced. Measurement of particle density in rapidity is a convenient way to describe 
heavy-ion collisions. It has also been suggested that fluctuations in pseudo-rapidity dis- 
tributions is a signature of phase transition from hadronic matter to Quark-Gluon plasma. 
Further, pseudo-rapidity density can be related to a thermodynamic quantity,entropy den- 
sity, in heavy-ion collisions. All these motivate us to study the pseudo-rapidity distribu- 
tions of photons ( dN-y/drf). We have studied dN 1 /drj in WA93 and WA98 experiments at 
CERN SPS []. These distributions for S + Au and Pb + Pb reactions for different centrality 
classes are shown in Figure ^. They are found to be good gaussians. The results from 
VENUS are also shown in the figures. One observes that VENUS under predicts data for 
central collisions. 

The best way to study the pseudo-rapidity distributions is to look at the shape parameters 
of the distributions. The shape parameters are the pseudo-rapidity density at mid-rapidity 
iPmax), width of the pseudo-rapidity distribution (er) and pseudo-rapidity peak {rj pea k). 
These are shown in Figure |^ as a function of transverse energy and number of partici- 
pating nucleons for WA93 and WA98 experiments respectively. The results have been 
compared to those obtained from VENUS event generator. The pseudo-rapidity density at 
mid-rapidity is found to increase with increase in centrality, which can be understood from 
simple geometrical picture of the collision. Results from VENUS is also found to follow 
the similar trend. For higher centrality of the reaction VENUS under predicts the data. 
The width of the pseudo-rapidity distributions for various centrality classes are found to be 
similar within the quoted errors and the trend is well explained by results from VENUS. 
The pseudo-rapidity peak for the WA98 experiment was found to be 2.92. 
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Figure 2. Pseudo-rapidity distributions of photons for S+Au reaction at 200-AGeV 
(left plot) and Pb+Pb reaction at 158-AGeV (right plot). The solid histograms are the 
corresponding distributions obtained from the VENUS event generator. 




"pr-t-T 



50 



100 150 200 

E T (GeV) 



250 



2 300 



200 
100 


1.8 
1.6 
1.4 
1.2 

1 
0.8 
0.6 



• Pb+Pb 
■ Pb+Nb 
▲ Pb+Ni 
Venus(Pb+Pb; 




50 100 150 200. .250 300 350 400 450 



part 



Figure 3. Pseudo-rapidity density (p ma x), and width of the pseudo-rapidity distri- 
bution (a) are shown as a function of centrality of the reaction (defined through the 
Et values for WA93 results (left plot) and number of participating nucleons for WA98 
results (right plot)). 
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3.2 Scaling of photons 

It is important to study the scaling of particle multiplicity as they test the various models 
for particle production. Also various experimental signatures require comparison of ob- 
servables of different system sizes, hence a proper understanding of scaling is essential. 
While scaling with number of collisions arises naturally in a picture of a superposition 
of nucleon-nucleon collisions, with possible modifications by initial state effects, the par- 
ticipant scaling is more naturally related to a system with strong final state re-scattering, 
where the incoming particles lose their memory and every participant contributes a similar 
amount of energy. The scaling behavior can therefore carry important information on the 
reaction dynamics. It is therefore of interest to study the scaling properties with respect 
to number of participants or collisions. We have studied the scaling of total number of 
photons with number of participating nucleons [] . The results for photons from Pb+Pb 
collisions at 158-AGeV are shown in Figure Fitting the data points to the function 
C x Np art , yields the value of a to be 1.12 ± 0.03. Similar analysis for charged particles 
in the same experiment and using the data from the silicon pad multiplicity detector gives a 
value of a = 1.07 ± 0.05 []. Within the quoted systematic errors the value of a are similar 
for both photons and charged particles and indicates a deviation from the picture of a naive 
Wounded Nucleon Model (a = 1). 




Figure 4. Scaling behavior of photons. Integrated number of photons are plotted as a 
function of number of participating nucleons. The solid lines show power-law fit to the 
data, which yields the value of exponent, a= 1.12 db 0.03. 



3.3 Event-by-event mean transverse momentum of photons 

Variation of < pr > of produced particles with global multiplicity could provide signature 
of phase transition from hadronic matter to Quark-Gluon plasma. The expected behavior 
will be similar to that of variation of entropy density with temperature, where temperature 
increases with increase in entropy density before phase transition then remains unchanged 
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Figure 5. The mean transverse momentum, < pr >, of photons as a function of 
centrality. For the S+Au collisions (left plot) centrality is defined by Ezdc / Eb ea m 
for the Pb+Pb collisions it is defined through the pseudo-rapidity density of photons at 
mid-rapidity,p maa; . The < pr > values obtained from the VENUS event generator are 
superimposed for comparison. 



during the transition and finally again increases with entropy density after phase transition. 
Entropy density is related to pseudo-rapidity density and transverse momentum to temper- 
ature. Hence a study of variation of transverse momentum with pseudo-rapidity density 
(or centrality of reaction) can be an experimental signature of QGP formation. 

We have measured event-by-event < pt > of photons in both WA93 [] and WA98 exper- 
iments [], using the relation (pr) = E T m / N 1 . Where E T rn is the transverse component of 
the electromagnetic energy measured in the mid-rapidity calorimeter, and N 7 is the num- 
ber of photons measured by PMD in a given common r)-<p region of the two detectors. The 
results are shown in the Figure ||. One observes that for both WA93 and WA98 experiment, 
the < pt > values are constant as a function of centrality of the reaction within the quoted 
systematic errors and agrees well with those obtained from VENUS. 



3.4 Event-by-event fluctuation in photon multiplicity 



Event-by-event fluctuations in global observables such as photon multiplicity has gener- 
ated lot of interest in recent years []. This is because fluctuations in the global observables 
are related to thermodynamical properties of matter, more specifically multiplicity fluctu- 
ations can be related to matter compressibility. Further studying the centrality and rapidity 
acceptance of fluctuations can help in finding the possible existence of the tri — critical 
point in the QCD phase diagram []. We have studied the centrality dependence and rapidity 
dependence of fluctuations in photons and the results are shown in Figure ^ []. The fluctu- 
ation in photon multiplicity (w 7 ) is defined by the relation u 7 = er 7 2 / < A^ 7 >. Where a 1 
is the standard deviation and < iV 7 > the mean of the gaussian photon multiplicity distri- 
bution for a particular centrality class. For results on the centrality dependence of photon 
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Figure 6. Left plot : Relative fluctuations, cj 7 of photons as a function of number 
of participants. These are compared to calculations from a participant model and those 
from VENUS event generator. Right plot : Photon multiplicity fluctuations for two r\ 
acceptance selections. The open squares represent estimated values of fluctuations in 
0.5 unit of Ar\ from the observed fluctuations in 1.0 unit A??. / is the ratio of mean 
photon multiplicity for 0.5 unit of Arj to that for 1.0 unit of Arj. 

multiplicity fluctuation the data has been compared to the results obtained from VENUS 
event generator and a simple participant model. The results of fluctuation in photon data 
are presented as a function of mean number of participant nucleon and are found to be in 
agreement with those form the model calculations within the quoted systematic errors. For 
the rapidity acceptance studies, the decrease in fluctuation with decrease in acceptance has 
been explained using a simple statistical model. The model assumes that the distribution 
of particles in a smaller rapidity window Ai] follows a binomial sampling []. The results 
have been presented as a function of 2% bins in cross section of total transverse energy. 

3.5 Formation of Disoriented chiral condensates at SPS 

The formation of hot and dense matter in high energy heavy-ion collisions has the possi- 
bility of creating a chiral symmetry restored phase in the laboratory. A consequence of this 
is the possibility of the formation of disoriented chiral condensate (DCC) []. The detec- 
tion and study of DCC is expected to provide valuable information about the chiral phase 
transition and vacuum structure of strong interactions. 

DCC formation in a given domain would be associated with large event-by-event fluc- 
tuation in the ratio of neutral to charged pions in a certain domain. Since neutral pions 
decays to photons and majority of the charged particles are charged pions. The experimen- 
tal signature of formation of DCC would be large event-by-event localised fluctuation in 
ratio of photons to charged particles. Following this theoretical predictions, WA98 recently 
carried out a detailed study of central Pb+Pb events and has been published in Ref. []. 
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laDie i. 


Type of fluctuations preserved by 


various 


mixed events 




Fluctuation 




Mixed 


Event 






Ml 


M2 


M3-ch 


M3-7 


7V 7 -only 


No 


Yes 


Yes 


No 


Nch ~ only 


No 


Yes 


No 


Yes 


correlated JV 7 -iV c h 


No 


No 


No 


No 



The data presented here were taken during December 1996 at CERN SPS with 
158 -A GeV Pb ion beam on a Pb target of thickness 213 /im. The Goliath magnet was 
kept off during these runs. A total of 85K central events, corresponding to the top 5% of 
the minimum bias cross section as determined from the total transverse energy measured 
by the mid-rapidity calorimeter (MIRAC), are analyzed. The analyzed data corresponds to 
a common -q (2.9 < i] < 3.75) and <f> coverage of charged particle and photon multiplicity 
detectors and for different segments in azimuth (<f>). 

The measured results are interpreted by comparison with simulated events and with sev- 
eral types of mixed events. Simulated events were generated using the VENUS 4.12 event 
generator with default parameters. The output was processed through a WA98 detector 
simulation package in the GEANT 3.21 framework. Due to the inherent uncertainties in 
the description of "normal" physics and detector response in the VENUS+GEANT sim- 
ulations, the observation of an experimental result which differs from the case with zero 
DCC fraction in simulation cannot be taken alone as evidence of DCC observation. For this 
reason four different types of mixed events have been created from the real events in order 
to search for non-statistical fluctuations by removing various correlations in a controlled 
manner while preserving the characteristics of the measured distributions as accurately as 
possible. The first type of mixed events (Ml), are generated by mixing hits in both the 
PMD and SPMD separately, with no two hits taken from the same event. Hits within a 
detector in the mixed events are not allowed to lie within the two track resolution of that 
detector. The second kind of mixed events (M2) are generated by mixing the unaltered 
PMD hits of one event with the unaltered SPMD hits of a different event. Intermediate be- 
tween the Ml and M2 kinds of mixed events is the case where the hits within the PMD are 
unaltered while the SPMD hits are mixed (M3-7), or the SPMD hits are unaltered while the 
PMD hits are mixed (M3-ch). In each type of mixed event the global (bin 1) -/V 7 _iik -A r c h 
correlation is maintained as in the real event. The specific physics issues probed by each 
mixed event are tabulated in Table 2. 

In order to look for any possible localized fluctuation in photon and charged particle 
multiplicities, which may have non-statistical origin, it is necessary to look at the cor- 
relation between iV 7 _iik c and at various segmentation in r\ and <f>. Event-by-event 
correlation between A 7 _iik c and A c h has been studied in 0-segments, by dividing the en- 
tire 0-space into 2, 4, 8 and 16 bins. Fig. [7] shows the scatter plots of A 7 _iik e and N c h 
distributions. The correlation plot for each of the segments, starting with the case of 1 
bin which corresponds to no segmentation are also shown. A common correlation axis 
(Z) has been obtained for the full distribution by fitting the A 7 _iik c and AT c j, correlation 
with a second order polynomial. The correlation axis with fit parameters are shown in the 
figures. The distances of separation {Dz) between the data points and the correlation axis 
have been calculated with the convention that Dz is positive for points below the Z-axis. 
The distribution of Dz represents the relative fluctuations of A 7 _i;k c and A c h from the 
correlation axis at any given <f> bin. In order to compare the fluctuations for different <j> bins 
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500 




Figure 7. Left plot : Scatter plot showing the correlation of N 1 _ likc and 7V c h for 
1,2,4,8 and 16 bins in azimuthal angle, <f>, for the common coverage of the PMD and 
SPMD. The common correlation axis (Z) obtained by fitting the distributions to a sec- 
ond order polynomial is also shown along with the fit parameters. Right plot : The Sz 
distributions for data (solid circles), Ml mixed events (solid histogram) and simulated 
events (dashed histogram) for various bins in azimuthal angle. Distributions of other 
mixed events are not shown for clarity of presentation. 



in the same level, we use a scaled variable, Sz = Dz / 's(Dz), where s(Dz) represents 
the rms deviation of the Dz distribution for VENUS events, obtained in a similar manner. 
The presence of events with localized fluctuations in jV 7 _iik and iV c h, at a given bin, is 
expected to result in a broader distribution of Sz compared to those for normal events at 
that particular bin. For comparison, Sz distributions are also obtained for mixed events. 

The Sz distributions calculated at different <f> bins are shown in Fig. ^ for data, Ml, 
and VENUS events. The distributions for other mixed events are not shown for clarity of 
presentation. The small differences in the Sz distributions have been quantified in terms of 
the corresponding rms deviations of these distributions as shown in Fig. ^| The statistical 
errors on the values are small and are within the sizes of the symbol. The bars represent 
statistical and systematic errors added in quadrature. The rms deviations of M2-type of 
mixed events are found to agree with those of the experimental data within errors. Thereby 
suggesting the absence of event-by-event localized correlated fluctuations in 7V 7 _iike and 
AT c h or DCC-type fluctuations. The rms deviations of Ml-type of mixed events are found 
to be lower than those obtained for data for 2, 4 and 8 bins in <f>. The results form M3-type 
of mixed events are found in between those obtained form data and Ml-type mixed events. 

These results indicate the presence of localized fluctuations in data, due to both photons 
and charged particles for certain bins in <fr. Within the context of a simple DCC model [], 
upper limits on the presence of localized non-statistical DCC-like fluctuations has been set 
to be 10" 2 for DCC domain size with Ac/) between 45-90° and 3 x 1(T 3 for A</> between 
90-135° with 90% confidence level. 
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Figure 8. The root mean square (rms) deviations of the Sz distributions for data, 
mixed events and simulated events as a function of number of bins in azimuthal angle. 



4. Photon multiplicity detectors at RHIC and LHC 

Photon multiplicity detectors will be installed in STAR experiment at RHIC and in ALICE 
at LHC. The detector technology is different from those used at SPS. PMD in STAR and 
ALICE will be a gas based detector with honeycomb cells []. Typical detector parameters 
for the PMD to be used in STAR and ALICE is given in the Table 3. 



Table 3. PMD in STAR and ALICE 



Basic features 



STAR-PMD 



ALICE-PMD 



Expected year of Installation 

Beam and target 

CMS Energy 

No. channels 

Readout 

Distance from target 
r) coverage 
Efficiency (central) 
Purity (central) 



2002 
Au + Au 
200 A GeV 
82944 

Gassiplex+CRAMS 

5.5 m 

2.3-3.9 

62% 

61% 



2005 - 2006 
Pb + Pb 
5.5 A TeV 
270,000 

Gassiplex+CRAMS 

3.5 m 

2.3-3.5 

64 % 

60% 



PMD will address all the physics issues discussed in Section 3. It will have overlap of 
about one unit in r) with full <p with the forward time projection chamber (FTPC) in STAR. 
The momentum information of charged particles in FTPC will strengthen the DCC-type 
of study at RHIC. This is demonstrated by a simple calculation given below. In addition, 
PMD will be able to give estimates of transverse electromagnetic energy in ALICE. High 
multiplicity environment at ALICE will be suitable for such estimates as the statistical 
fluctuations on transverse electromagnetic energy measurements gets reduced. 
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4.1 DCC search in STAR 

The pr information of charged particles will enhance the possibility of DCC search and 
verify the various features of DCC domain, such as DCC-pions are low px pions and DCC 
formation may lead to low px enhancement in pions. The photon multiplicity detector in- 
conjunction with a charged particle detectors which gives multiplicity along with particle- 
wise momentum information, such as the FTPC in STAR experiment can be used. 

In order to show the utility of px information we carried out the following simulation 
study. The basic detector setup for DCC study was implemented by putting a photon mul- 
tiplicity detector and a charged particle multiplicity detector with momentum information, 
having a common coverage (as in STAR ) with realistic detector parameters. The photon 
counting efficiency was taken as 70% ± 5% and contamination in photon sample was taken 
to be 30% ± 5%. The detectors accepts particles with momentum above 30 MeV/c. The 
efficiency of charged particle detector was taken as 90% ± 5% with px resolution of =E^- = 
0.2. Then a DCC-type fluctuations in low px pions was introduced through a simple DCC 
model as described in Ref. []. Events with different fractions of DCC-type events were 
generated by mixing normal events and DCC-type events. Then the resultant events were 
analysed using DWT analysis. The analysis was done with and without px cut on charged 
particle data(pr i 150 MeV/c). The results are shown in Figure^. One finds the strength, 
(, of the DCC signal as a function of percentage of events being DCC-type is higher with 
px cut than without px cut. This clearly demonstrated that with px information of charged 
particles the sensitivity to search for DCC at RHIC increases. 
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0.2 r 
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Figure 9. Variation of £ as a function of % of events being DCC-type, for charged 
particle detector with pr information. Also shown are the corresponding results with 
charged particle detector without px information. 
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5. Summary 

A high granularity preshower photon multiplicity detector was successfully implemented 
in the WA93 and WA98 experiments at CERN SPS. The PMD provided the first ever 
measurement of multiplicity, rapidity and azimuthal distributions of photons in ultra- 
relativistic heavy-ion collisions. Analysis of these distributions showed that photon multi- 
plicity is an important global observable, providing significant information regarding parti- 
cle production and reaction mechanism in heavy-ion collisions at relativistic energies. The 
study of pseudo-rapidity distributions of photons, event-by-event measurement of mean 
transverse momentum of photons, scaling of photons and event-by-event fluctuations in 
photon multiplicity suggest that it is complementary to charged particle measurements. 
In addition to this, correlation between photon and charged particle multiplicity has been 
found to provide important information regarding formation of DCC in relativistic heavy- 
ion collisions. A gas based PMD will be implemented in the STAR experiment at RHIC 
and in ALICE experiment at LHC. It will address all the physics issues, which were studied 
using the PMD at CERN SPS. However, DCC search will get a boost due to momentum 
information of charged particles at STAR. With very high multiplicity environment at AL- 
ICE, PMD can provide estimates of transverse electromagnetic energy at LHC. 
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